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Abstract 

Both uniaxial and biaxial nematic liquid crystals are defined by orientational ordering of their 
building blocks. While uniaxial nematics only orient the long molecular axis, biaxial order implies 
local order along three axes. As the natural degree of biaxiality and the associated frame, that can 
be extracted from the tensorial description of the nematic order, vanishes in the uniaxial phase, we 
extend the nematic director to a full biaxial frame by making use of a singular value decomposition 
of the gradient of the director field instead. New defects and degrees of freedom are unveiled and 
the similarities and differences between the uniaxial and biaxial phase are analyzed by applying 
the algebraic rules of the quaternion group to the uniaxial phase. 



Nematic liquid crystals are revered for their geometrically complex and visually com- 
pelling defect structures, stabilized by topological constraints pp. The elementary rules of 
homotopy that govern these defects [2HS] imply ambiguities in defect classification when 
disclination lines are involved due to the action of the first homotopy group on itself and 
on the second homotopy group. Nematic defects have been probed in both nematic and 
cholesteric colloids and emulsions p4T2]. highly confined geometries [T314T5] . and optically 
manipulated liquid crystals [161 El- So robust are they, that they can be manipulated to 
reproducibly form linked or knotted disclination lines [TH]-EI] as well as other topologically 
interesting objects Recall that uniaxial nematics can be interpreted as a highly sym- 

metric special case of biaxial nematics [U [6] , suggesting an opportunity to study nematic 
defects with tools that are not available in the standard homotopy theory. In this article, we 
explore the similarity between distortion patterns found in uniaxial nematic fields and de- 
fects in biaxial phases by introducing a new biaxial frame derived entirely from deformations 
of the uniaxial director field n. Like the Frenet-Serret frame of a curve or the principal axes 
frame of a surface [21], our new frame has a well-defined (differential) geometric meaning 
and allows us to provide a topological characterization to the director geometry. The "quasi- 
defects" in this new frame allow us to apply the well-developed theory of biaxial nematics 
and to include the non-topological "escaped defect" [25] in our classification, embellishing, 
for instance, our understanding of the double-twist tube construction [26j of the blue phases. 
A motivation for this investigation is the study of similar quasi-defect structures in optics, 
where topological filaments in the derivative of a complex scalar field determine the topol- 
ogy of optical vortices [271 EH]- We demonstrate our technique on numerical models of blue 
phases and discuss the implications of newly extracted information. 

Uniaxial nematics consist of elongated non-polar molecules that tend to align in a par- 
ticular direction in space, taking directions in the manifold MP^. The director is specified 
by a unit vector n up to sign. As a result, uniaxial nematics accommodate both line defects 
(disclinations) and point defects [H E] since 7ri(]RP^) = Z2 and 7r2(]RP^) = Z. While point 
defects can be oriented, a line defect, a disclination line, winds the director by vr leading 
to a sign inconsistency. Recall, however, that non-defect states exist as disclinations with 
a non-topological integer winding number. The famous escape into the third dimension 
renders these smooth [2S] in n. As a line field alone, the nematic director does not have 
an intrinsic biaxial nature. In order to define a biaxial structure, we turn to the tensor of 
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gradients, diUj, which provides the additional structure necessary to define an entire frame. 
Multiplying from the left by a unit vector extracts a directional derivative of n; there are 
special orthogonal directions [w^,w^, w^], in which the magnitudes ||w7°c?jnj|p of the deriva- 
tives are extreme, which can also be seen as an eigenvalue problem. The derivatives in these 



a different orthonormal frame in the domain of the matrix in such a way, that the singular 
values (Jq, are positive semi-definite and satisfy o"i > o"2 > 0-3 > 0. Because the director n is 
a unit vector, it is in the kernel of the derivative tensor, {dinj)nj = 0, which implies = 
and = n. This makes the axes and orthogonal to n and decorate the uniaxial phase 
with a biaxial order! We thus decompose the gradient tensor as: 



By construction, the derivative of the director is largest in the direction w , changing towards 
with the rate of cxi. The derivative in the direction points towards with a lower rate 

of (72. The remaining vector w'^ = w marks the direction in which the director is constant. 

This is known as the singular value decomposition (SVD) of the matrix diTij and proves to 

be just the thing we need. 

It is amusing to note that the derivative tensor is intimately connected to the topological 

hedgehog charge through the Gauss integral: [HI [29] , 



Thus the streamlines of the vector field w = aia2W trace the preimage of a particular 
director orientation as a three-dimensional generalization of the schlieren texture [30j and 
the topological charge is simply the fiux of these streamlines through the enclosing surface. 
The streamlines can only terminate on the singular nematic defects or where the w field 
disappears, i.e., when (T2 = 0. The latter are saddles in the derivative field, where the 
director is constant in two directions. 

The SVD frame decorates the director field with two new vectors that encode the trans- 
verse degrees of freedom and can be interpreted as a coordinate in SO (3) / D2 to parameterize 
biaxial order. In regions where diUj does not vanish, the frame is continuous everywhere 



directions take the form of wfdiUj = aan'j for each a G {1,2,3}, where [n^,n^,n' 



.^] define 
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except on a one-dimensional set of points, which we collectively call line defects or disclina- 
tions in the SVD triad, as they are analogous to disclinations in biaxial nematics. There are 
two types of defects in this frame. On native disclinations - line and point defects in n (for 
simplicity, we can treat the point defects as small disclination loops [8J), the derivative tensor 
diverges and the frames and singular values are ill-defined. Away from the defects in n, the 
SVD decomposition can still give a frame with degenerate axes wherever two singular values 
coincide; we call these locations quasi- disclinations. Both the n-frame and the w-frame share 
the defects, which are characterized by singularities and degeneracies in the singular values. 
Each disclination line is characterized by a pair of degenerate axes and by the amount of 
rotation of these axes around the remaining nondegenerate axis - the winding number of the 
disclination. We will draw an analogy with the biaxial phase and its first homotopy group - 
the quaternion group written in terms of unit quaternions, {1,-1, ±i, ±j, ilk} |8], [31]. For 
the disclinations with half-integer winding number, we shall assign a unit quaternion i to a 
defect in axes {n^, n}, j to a defect in {n, n^} and Ik to a defect in {n^, n^}. All disclinations 
with odd integer winding number belong to the class —1, regardless of which pair of axes 
they involve. Finally, the disclinations with even integer winding number are all trivial, in 
the 1 class. 

First, we study the native nematic disclinations. They all involve the director n as one 
of the degenerate axes. Moreover, on a tight circle around the disclination, the director 
rotates very rapidly in the plane defined by the two degenerate axes, so the singular value 
associated with this direction diverges at the defect core. Since the singular values are sorted 
by magnitude, all native nematic disclinations are, by construction, degenerate in the axes 
{n, n^}. To see this, consider a plane perpendicular to the defect line. In this plane n winds 
around ever more rapidly as we approach the core. By definition, this rapid winding is into 
the direction at all points, wjdirij = crin^ with ai diverging. Note that since n ■ = 0, 
wldiU^ = —aiUj + 7n| where 7 is finite. It follows that the winding is between n and 
and that the defects associated with the degeneracy of the smallest two singular values, 
(J2 = 0, must be topologically trivial in the n frame. The disclinations with a half-integer 
winding number, already known from the homotopy theory of the bare director field, now 
also have the perpendicular axis performing a half- integer turn (Fig. [l^). The addition of 
the perpendicular axes also reveals disclinations with an odd integer winding number, which 
are not topologically distinguished in the standard uniaxial setting, but are well-defined in 
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FIG. 1: (color) Sketched profiles for both the native disclinations (a,b) and the hidden quasi- 
disclinations (c,d). The native nematic disclinations involve half-integer (a) or integer (b) rotation 
of the director n and the perpendicular axis n^. The quasi-disclinations are similar, except the 
auxiliary axes {n^,n'^} are degenerate instead of {n, n^}. Each depicted structure is just one of 
the many representatives of its class, as continuous rotations preserve the topology of the defects. 
(e,f) Numerical model of blue phases I and II with isosurfaces representing half-integer native (red) 
and quasi- (cyan) disclinations. The blue phase I has the quasi-disclinations running along the axes 
of double twist cylinders, extending infinitely along each lattice direction. The blue phase II has 
four native disclinations and four quasi-disclinations meeting at the same point, spanning all the 
diagonals of the cube in an alternating order, (g) At the quasi-disclinations, the axes (red) and 
(green) rotate by —n around n (blue), (h) At the native disclinations, the n and axes show 
a similar behavior. (f,g,h) In the blue phase II, exchanging the director n with converts one 
disclination into the other and effectively rotates the unit cell by 90 degrees. 

our definition (Fig. [T]d). 

The line-like quasi-disclinations, on the other hand, have smooth director complexions, 
and so the degenerate axes with nonzero winding number must be the invisible perpendicular 
axes {n^,n^}. These disclinations are located where the singular values (Xi = a2, precisely 
when we can no longer distinguish the two axes in the decomposition ([T]) . These defects can 
have either half-integer or integer winding number (Fig. [lj3,d). 

The odd integer- winding- number disclinations all belong to the same —1 class of biaxial 
disclinations and can transform one into the other: an unescaped integer winding number 



disclination with a singular core in n can escape into the third dimension, which removes the 
singularity, but still shows up as a defect in the perpendicular axes {n^,n^} of the frame. 
Our construction therefore unambiguously locates both the escaped and unescaped integer 
nematic disclinations, which would be impossible to locate from the director field alone. 

What about point defects in the uniaxial phase? Biaxial nematics cannot have point de- 
fects as 7r2[SO(3)/-D2] = 0. As we mentioned, to consider the biaxial structure we will inflate 
all point defects into small disclination loops carrying the nontrivial element of 7ri(]RP^). 
Each native disclination loop can either carry an odd or even topological point charge, mea- 
sured by the second homotopy group, and it can be linked by an even or an odd number 
of other disclination loops [3, El |32l |33] . The half-integer native disclinations and all flavors 
of integer disclinations form an abelian subgroup of the quaternion group j'^ with z/ G Z4. 
This periodicity of four, consistent with the theoretical result given by the torus homotopy 
group [U [32] is seen in the specialized form for —1/2 disclination loops and their self-linking 
numbers [31], and in the generalization to disclination loops with arbitrary cross section 
|35j . The orientation of the disclination profile in the immediate neighborhood of the sin- 
gularity holds the information about the two degenerate axes that vary rapidly around the 
singularity, foreshadowing the significance of quaternions and differential definition of our 
biaxial framing. The statement that the nematic hedgehog charge is a residual of linking a 
— 1 biaxial disclination |8] is intrinsic in our construction. 

Unlike in proper biaxial nematics, the point defects are still present in our system - 
as small native disclination loops, linked by quasi-disclinations. To further explore the 
connection between the point charge and the threading of the native loops, recall that 
the topological charge q in nematics is the degree of mapping — )■ MP^ from a closed 
measuring surface in the nematic medium to the unit sphere. The director n on this enclosing 
surface is in fact a g-covered sphere, which by the Riemann-Hurwitz theorem has the Euler 
characteristic of x = 2g. The transverse axes {n^,n^} form a tangent plane to MP^ and 
so, by the Poincare-Brouwer theorem, the total winding number of their surface defects 
must be equal to the Euler characteristic. Every sphere that encloses a point defect is 
penetrated by quasi-disclinations, with the sum of their winding numbers equal to 2q. All 
quasi-disclinations terminate on the point defects and thread native disclination loops that 
carry a topological charge. 

Similarly to the native disclinations, the half-integer quasi-disclinations form a distinct 
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k'^ subgroup, with the same algebraic structure and hnking rules that we found for the 
native nematic disclinations. However, the coexistence of half-integer disclinations of both 
flavors reveals the nonabelian nature of the quaternion group: the j and Ik disclination loops 
cannot link without creating a —1 line connecting the loops [36j. An additional feature of 
our gradient framing is the absence of i quaternions representing the impossible disclinations 
with a degenerate pair of axes {n^,n}, as discussed above. This is a variation of Poenaru's 
theorem [37] that limits the merging of defects in gradient fields, as in smectics [38]. Were 
a half-integer quasi-disclination (k) to merge with a half- integer native disclination (j), 
the resulting defect would have the impossible signature i, so the quasi-and native defects 
naturally avoid each other. They can only meet at discrete points - at point defects, which 
we have shown act as sources or sinks for the quasi-disclinations. As in systems with broken 
translational symmetry [H |38] , the missing quaternion causes the homotopic description to 
be incomplete in this case: the ±i elements of the fundamental group have no realization in 
the sample. 

Blue phases are a suitable system for studying the biaxial defects on practical data, 
as they show nonuniform behavior without complicated boundary conditions. We use the 
finite difference relaxation method, based on the Landau-de Gennes model used in Ref . [39] . 
The Q-tensor field was used to retrieve the order parameter S, and the director field n, 
which was subsequently differentiated and decomposed with SVD, giving singular values 
and framing information for each point in space. A large resolution of 80 points along each 
direction was used, as the positions of the singularities are sensitive to errors caused by finite 
difference approximation of the derivatives. The zeroes of the order parameter 5* were used 
to determine the position of the native disclinations, while the zeroes of a functional ai — 02 
were used to find the quasi-disclinations. 

Both blue phases form a periodic cubic lattice. The blue phase I consists of straight 
native disclinations, extending infinitely in the direction of body diagonals and offset by 
half of the cell spacing to avoid each other. The rest of the bulk can be roughly explained 
as three mutually perpendicular double-twist cylinders, extending in the direction of main 
coordinate axes |3lJ. Plotting the near-zero isosurfaces of o\ — 02 reveals three mutually 
perpendicular infinitely extending quasi-disclinations that approximately follow the double- 
twist cylinders, which is not unexpected, as a\ = 02 condition implies the rate of change 
is equal in two directions, which is likely to occur near the axis of a double twist cylinder 
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(Fig. [Tp). In particular, the geometric pattern of the cross section has a half-integer winding 
number and looks like a characteristic three- fold profile of a —1/2 nematic disclination line 
in the axes and n^. 

The blue phase II consists of a cross-linked network of disclinations. The unit cell contains 
two junctions where four diagonal native disclination lines meet in a tetrahedral formation. 
The quasi-disclinations also extend from one junction to the other in straight diagonal lines, 
forming a tetrahedral structure, dual to the native one (Fig. [T]F). The junctions are thus 
highly degenerate, with four native and four half-integer quasi-disclinations meeting there 
and extending to all 8 vertices of the unit cube. An inspection of the disclination cross 
sections again reveals that both the native and the quasi-disclinations have a three-fold 
profile (Fig. [l^,h). Furthermore, the cross sections of both disclination types are exact 
copies of each other, with the axes n and in exchanged roles. In fact, the framing 
across the entire space possesses such a symmetry, that an exchange of these axes has the 
same effect as a rotation of the unit cell by 90°, which implies that even though the point 
symmetry group around the central point is the tetrahedral group, there are hints of the full 
cubic symmetry in the system. The axis in fact has a full cubic symmetry, up to a small 
perturbation that depends on the difference in the free energy costs of native and quasi- 
defects and could therefore be used as a model director field for a periodic cubic structure 
with an eight-way central junction of half- integer disclinations [10] . In general, every director 
field has adjoint fields and with equal or higher symmetry as the original director, 
which have a potential use as model director fields for related problems, as they already 
respect the boundary conditions and approximate energy constraints. 

In this paper, we explored the uniaxial nematic as a special kind of biaxial nematic with 
hidden perpendicular degrees of freedom. Instead of using the eigenvectors of the tensorial 
order parameter, which is highly degenerate in the uniaxial phase, we retrieved the missing 
perpendicular axes by using the SVD decomposition to extract a smoothly varying frame 
from the director derivative, treating the disclination lines as simple defects in the biaxial 
frame. This frame is nevertheless related to the frame, retrieved from the Q-tensor, as the 
first encodes the spatial variations of the director and the latter describes thermal fluctua- 
tions of the molecular director. Continuity of the director field and integrability conditions 
near singular defects give rise to rules that restrict the set of allowed defects, resulting in 
an intricate structure that is similar, yet not equivalent to that of general biaxial defects. 
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Beside the native uniaxial nematic disclinations, known from the conventional homotopy 
analysis, we uncover another type of disclinations that arise through hidden biaxial order - 
topologically unavoidable patterns in the elastic response, such as the escaped integer discli- 
nation lines, which were impossible to treat topologically in bare nematic director fields, but 
can now be computationally detected. With the elusive defects pinpointed deterministically, 
the complete set of linking rules for disclination lines emerges and unifies the line and point 
defects under the same formalism. 

The SVD decomposition has a potential use in numerical and analytic calculations, as it 
connects the free energy with topology via the singular values. The n-frame we focused on in 
this paper, can be taken as a convenient choice of frame in the Mermin-Ho construction |2l] . 
Additionally, the conjugate w-frame also constitutes a biaxial frame that can be investigated 
in the future. As a visualization technique, the singular values help to locate the escaped 
defects and other characteristic features without resorting to visualization of vector fields. 

The technique illuminates an intricate link between geometry and topology independently 
of the physical meaning of the order parameter. We can find an underlying biaxial field and 
extract hidden degrees of freedom for a wide variety of materials that allow parametrization 
by an orthogonal frame. Of particular interest for future research are patterns and defects 
in chiral nematics, smectic liquid crystals and fields of polarized light, building on the 
connection between local differential structure and the traditional homotopy of defects. 
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